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determined for analytic surfaces and reasonably accurate results 
(maximum error of less than 20 percent) for ratios of the order 
of twice this limit. This conclusion is significant because some 
of the previous investigators [7], [10] have suggested that the 
method is not applicable to such nonanalytic surfaces. 
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Calculation of the Radiation Resistance of Loop Antennas 
with Sinusoidal Current Distribution 


A. RICHTSCHEID 


Abstvact—Loop antennas with dimensions comparable to the wave- 
length are sometimes of practical interest. In order to use such antennas, 
design engineers would need simple formulas or data about the radiation 
resistance of such antennas. For the reason, that nothing else has been 
known, the radiation resistance has been calculated by means of a digital 
computer, 


I. INTRODUCTION 


In normal radio communication, loop antennas are seldom 
used, and in most of these cases they serve only for receiving 
cases. The reasons are that mostly loop antennas have very little 
efficiency, especially when they are driven with uniform current 
(i.e., the circumference of the loop is small compared with the 
wavelength). It would be better, especially for all transmitting 
cases, if the circumference of the loop were comparable with the 
wavelength. The radiation resistance of the loop would be much 
greater than the ohmic losses and, therefore, the efficiency would 
increase. Unfortunately a solution of the radiation resistance of 
loop antennas with nonuniform current is not known. Therefore 
the radiation resistance of loop antennas with sinusoidal current 
distribution has been calculated (sinusoidal current distribution is 
a reasonable assumption for loops with a circumference up to 
about one wavelength). 


If. RapIATION RESISTANCE OF Loop ANTENNAS WITH UNIFORM 
CURRENT 


The radiation resistance of loop antennas with uniform current 
and dimensions that are small compared to the wavelength can 
easily be calculated and is well known [1]. The solution for this 


Manuscript received September J, 1975; revised April 26, 1976. | 
The author is with the Fachgebiet Hochfrequenztechnik, Technische 
Hochschule Darmstadt, D 6100 Darmstadt, West Germany. 


889 


[8] S. R. Murphy and G. E. Lord, “Scattering from a sinusoidal surface 
—A direct comparison of the results of March and Uretsky,” J. 
Acoust. Soc. Amer., vol. 36, pp. 1598-1599, 1964. 

[9] R. Petit and M. Cadilhac, ‘Sur la diffraction d’une onde plane par 
un réseau infiniment Son gucieHe. C. R. Acad. Sci. (Paris), ser. A-B, 
vol. 262, pp. 468-471, 966. 

{10] R. F. Millar, “On the Ravicigh assumption in scattering by a periodic 
surface II,” Proc. Camb. Phil. Soc., vol. 69, pp. 217-225, 1970. 

[11] G.R. Jiracek, * ‘Numerical comparisons of a modified a roneh with 
other rough surface EM scattering solutions,” JEEE Trans. Antennas 
Propagat. (Commun.), vol. AP-21, pp. 393-396, May 1973. 

{12] H. Ikuno and K. Yasuura, “Improved point- -matching method with 


application to scattering from a periodic surface,” JEEE Trans. 
Antennas Propagat., vol. AP-21, pp. 657-662, Sept. 1973. 
[13] s “Numerical analysis of diffraction gratings,” Trans. [ECE 


Japan, vol. 54-B, pp. 182-189, 1971. 

[14] V. Jamnejad-Dailami, R. Mittra, and T. Itoh, “A comparative study 
of the Rayleigh hypothesis and analytic continuation methods as 
applied to sinusoidal grating,” JEEE Trans. Antennas Propagat. 
(Commun.), vol. AP-20, pp. 392-394, May 1972. 

[15] A. Wirgin and R. Deleuil, “Theoretical and experimental investigation 
ae be ed of blazed grating,” J. Opt. Soc. Am., vol. 59, p. 1348, 


case is 


= 197 ()* (1) 
A 


where 


R, radiation resistance 
/ circumference of the loop 
A wavelength. 


Equation (1) is valid for all cases where //2 < 1/4, i.e., where the 
current distribution is not differing essentially from the constant 
case. 

In [2] the radiation resistance of loop antennas with uniform 
current and, compared to the wavelength, great dimensions, has 
been calculated. But this does not solve the problem; because in 
real physics, there will necessarily result an—theoretically—ideal 
sinusoidal current distribution. Sometimes there is proposed, to 
enforce a uniform current distribution by means of multiple 
feeding. In most cases this is not practicable since there would be 
paralleled much little antenna sections. Each of these sections 
has by reason of the dividing, a radiation resistance, which is 
much smaller than the radiation resistance of the whole antenna. 
Therefore the terminal resistance will become very small, and it 
will be difficult to match the antenna. Moreover the efficiency 
will decrease by reasons of the long and numerous section 
feeders. In Fig. 1 there is shown the calculated radiation resistance 
after [2] in the case of a uniform current distribution. 


III. SrnusoIpDAL CURRENT DISTRIBUTION 


The radiation resistance R, of a thin ideal conducting loop 
antenna can be determined with the radiated power N and the 
current amplitude Jy by 


N = 44g? R,. (2) 
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Fig. 1. 


Radiation resistance of loop antenna with uniform current 
distribution. 


current element 


Fig. 2. Orientation of loop antenna. 


. The true power density P, which is radiated in radial direction, 
can be determined with Poyntings vector by 


1 
P=%3Ex H* =-— Ex E* 3 
z 5 GB) 


ce) 


E - electric field intensity 

ff magnetic field intensity 

H* conjugate complex value of H 

Zo 1202, characteristic impedance of the free space. 


For a more simple calculation it is advantageous to substitute 
.the loop antenna by a chain of many little elementary dipoles 
which have all constant, but differing from the neighbor current 
distribution. By summing up the field intensities of these elemen- 
tary dipoles one gets the field distribution of the whole loop 
antenna. Integrating Poyntings vector over a spherical surface 
will deliver the radiated power. 

The analysis procedure for getting the final integral whose 
solution delivers the radiation resistance is well known and often 
described in antenna literature. Therefore, the procedure of 
getting the integral is only shortly referred. The loop antenna 
‘shall be orientated as it is shown in Fig. 2. 
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The current distribution is described by 


Io cos (v i) 


Ip cos (2= ‘ v) (4) 


Tw) 


he 


where a is the radius of the loop and / is the circumference of the 
loop. (For loops with dimensions much greater than one wave- 
length this assumption is quite theoretical, nevertheless the 
results were calculated till //A = 10 and are plotted till //A = 2.5. 
This may better help to show the differences between the two 
modes). 

The loop is divided in many elementary dipoles with the length 
di, and the resulting field components E, and E, are determined. 
One elementary dipole is marked in Fig. 2. 

Considering the variation of the angle yw and the varying 
distances r of the current elements from the regarded coordinates 
of the point P (ro, 4,8), one gets for the far-field from one element 


Iga cos (lw/a) 


ie a Z sin (¢ — ¥) cos 6 
. a gika sin 6 cos (o-W) dw (5) 
ro 
_ io Wh]h 
dE, = —i2, 708 cos (¢ — ¥) 
e Jkro 


elka sin 6 cos (@—p) dw. (6) 
Fo 


To get the complete field, it is necessary to integrate all particular 
field-components. 
The power density P is delivered by 
1 
pone E,E,* + : 
2 2Z, 


0 a 


E,E4". (7) 


With (7) one can destinate whole the power N which is radiated 


bpd 


ra /2 
= ed [ {; {E,E,* + E,E,* hr? sin 8 de dg. (8) 
Zo g=0 Jb=0 
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With (2)-(7) one gets, after some conversions, 


2 ria ef 
R= Z, (;) { [ (Se 4 Ie tS FY 
4} Je=0 Je=0 
-sin@d0d¢ (9) 


wherein J; are integrals of the form 


Ps | : 
: _ - cos (ka sin 8 cos (¢@ — W)) dy 
ie f cos (kay) sin (¢ — W) cos 8 
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-sin (ka sin 8 cos (¢ — W)) dw. (10) 
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Fig. 3. Radiation resistance of loop antenna with sinusoidal current 


distribution versus constant current distribution. 


The solution of this triple integral is leading already after the 
first integral to Bessel functions of the order L£a+ jp), Where a 
and b may not be integer. A retrospective view at the way of 
solution in [2] shows, that the author has calculated just the 
last problem, which can be solved closed. 

In [2] it has been possible, to convert the Bessel functions (for 
constant current their order is integer) into sums, and therefore 
they remain integrable, not so in this special case. For sinusoidal 
current we have solved the integral (9) with a digital computer. 


IV. CALCULATION OF THE RADIATION RESISTANCE 
UsING A COMPUTER 


By solving the problem with a digital computer special atten- 
tion had to be paid in order to get sufficient accuracy. Normally 
the accuracy of the result will increase with an increasing number 
of probes. On the other hand, the accuracy will decrease if there 
are too many probes because of the accumulation of rounding 
errors and other inaccuracies of the computer. An optimal num- 
ber of probes had to be found. It has been shown, that the ac- 
curacy is strongly dependent on the number of probes in yw. Less 
critical are the other dimensions ¢ and @. Moreover it has to be 
observed, that there will be no interfering between probe number 
and wave number which would cause incorrect extinction of some 
field components. Therefore the probe distances have been 
chosen not constant, but increasing with a special function. The 
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resulting faults can be eliminated by applying a corregating 
weighting function. 

In.Fig. 1 the calculated radiation resistance is shown versus the 
wave number around the loop (dotted line). If one compares this 
with the other curve, it can be seen, that the radiation resistance 
for sinusoidal current distribution is expectively much less. In 
the region where it is //2 « 1, both curves are approximately 
equal. There also appears an oscillation around a medium value 
(for //A > 1 compare in [2]), but now this oscillation has a period 
of 1 in contrast to the period z. This leads to the explanation 
that in the case of a Joop antenna with uniform current distri- 
bution the oscillation is the effect of adding or subtracting field 
components, which are excited from diametric antenna segments. 
In the case of sinusoidal current distribution the oscillating is 
caused by adding or subtracting field components, which now 
result from neighboring segments. This effect is well known from 
linear antennas. 


YV. CONCLUSION 


It has been shown, that the calculation of the radiation resist- 
ance of a loop antenna with sinusoidal current distribution 
cannot be solved closed. Therefore the solution of the problem 
has been carried out by the aid of a digital computer. 
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Measured Field of a Directional Antenna Submerged in a Lake 
L. C. SHEN, R. W. P. KING, AND R. M. SORBELLO 


Abstract—Field patterns of a VHF directional antenna submerged in a 
lake are shown with the receiver both under and over the surface of the 
water, The measured relative strength of the electric field as a function 
of the radial distance from the source is also shown. The rate of decrease 
with distance of the radial component of the field in the water agrees with 
theoretical predictions. 


I. INTRODUCTION 


The transmission of electromagnetic waves over the surface of 
the earth has been a subject of scientific interest and technical 
importance since the formulation of the basic theory by Sommer- 
feld in 1909. At that time and for many years the primary 
motivation was radio communication between transmitting and 
receiving antennas erected above the earth. The application of 
the theory to subsurface communication between two sub- 
marines began with a thesis by Moore in 1951 and two un- 
published reports by Bafios and Wesley in 1953 and 1954. 
Expanded and corrected treatments are in the more recent 
literature [1]-[3]. In particular, the comprehensive work of 
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